Draft version November 12, 2012 

Preprint typeset using L^T^X style cmulateapj v. 5/2/11 



SHAPE EVOLUTION OF MASSIVE EARLY-TYPE GALAXIES: CONFIRMATION OF INCREASED DISK 

PREVALENCE AT Z > 1 

Yu-Yen Chang 1 , Arjen van der Wel 1 , Hans- Walter Rix 1 , Stun Wuyts 2 , Stefano Zibetti 3 , Balasubramanian 

Ramkumar 1 , and Bradford Holden 4 

Draft version November 12, 2012 



ABSTRACT 

We use high-resolution if-band VLT/HAWK-I imaging over 0.25 square degrees to study the struc- 
tural evolution of massive early-type galaxies since z ~ 2. Mass-selected samples, complete down to 
log(M I 'Mq) ~ 10.7 such that 'typical' (L*) galaxies are included at all redshifts, are drawn from pre- 
existing photometric redshift surveys. We then separated the samples into different redshift slices and 
classify them as late- or early-type galaxies on the basis of their specific star-formation rate. Axis-ratio 
measurements for the ^400 early- type galaxies in the redshift range 0.6 < z < 1.8 are accurate to 
0.1 or better. The projected axis-ratio distributions are then compared with lower redshift samples. 
We find strong evidence for evolution of the population properties: early-type galaxies at z > 1 are, 
on average, flatter than at z < 1 and the median projected axis ratio at a fixed mass decreases with 
redshift. However, we also find that at all epochs z < 2 the very most massive early-type galaxies 
(log(M /Mq) > 11.3) are the roundest, with a pronounced lack among them of galaxies that are flat 
in projection. Merging is a plausible mechanism that can explain both results: at all epochs merging 
is required for early-type galaxies to grow beyond log{M / Mq) ~ 11.3, and all early types over time 
gradually and partially loose their disk-like characteristics. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: structure — galaxies: ellip- 
tical and lenticular, cD — cosmology: observations 



1. INTRODUCTION 

In our theory of galaxy formation, the initial angu- 
lar momentum, radiative energy loss of the gas, and 
some degree of angular momentum conservation result 
in gas settling into disks before most of the stars form, 
which makes for oblate, rotating stellar systems. Ob- 
servationally, all but the most massive galaxies in the 
present-d ay universe have a disk-like structure and are 
rotating (iDresslerl 119801: iKormendv fc Di or govskil 119891: 
Marinoni et al.lll999l:lKrainovic et alll2008HHolden et ail 
2009; E msellem et alJl20lH iCheng et alll201lD . 

Nonetheless, even in disk-like, rotating galaxies, many 
stars reside in pressure-supported bulges. Photo- 
metric decompositions indicate that in the present- 
day Universe 58±7% of stars are in the sp heroids 
and 42±7% are in the disks ([Benson et all 120071 ). 
Violent mergers are thought to scramble the or- 
bits that orig i nally were formed a n d lived in disks 
(iToomrd [19771: IWhite fe Reesl [l97l IJoseph fc Wrightl 



iBernardi et~aIll20lH lEmsellem et aT1l2011t Hidden et al.l 



19851 ISchweizer et al l 119901: iKauffmann et al.l 119931: 
Steinmetz fc Navarro! [2002). This process dominates 
the evolution of the most massive galaxies, which 
do not show evidence for disks and are generally 
round an d entirely supported by p ressure instead of 
rotat io n (IJorgensen fc Franxl 119941 : IVincent fc Rvdeiil 
[20051 Ivan der Wel et all ([200911 hereafter vdW09); 
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([2012J, hereafter H12)). 

In this paper we use a low specific star formation as the 
definition of early-type galaxy, motivated by the smooth 
appearance of the light profiles of galaxies with little or 
no young stars or (star-forming) gas. Early-type galax- 
ies show a relatively abrupt change in their structure 
as a function of galaxy mass. vdW09 and H12 show 
that galaxies with log(M/M Q ) < 11 have a broad pro- 
jected axis-ratio distribution, indicative of a disk-like 
stellar body (with typical short-to-long intrinsic axis ra- 
tio about 1:3). In contrast, early-type galaxies with 
log(M/M Q ) > 11 rarely appear highly flattened, sug- 
gesting that their formation channel destroyed any pre- 
existing stellar disks, with 'dry' merging between gener- 
ally gas-poor progenitors as a very plausible mechanism. 

Such 'dry' mergers can reconcile the relatively late, and 
continuous, assembly of massive early types with their 
old, passively evolving stellar populations. In addition, 
the observed small radii and hig h densities of early-type 
galaxies at redshifts z > 1 (e.g., [van Dokkum et aD[2008l : 
van der Wel et all [2 0081 can also be explained by 'dry' 



merging (e.g. 



iKhochfar fc SiUd l2006t Ivan der Wel"et~atl 
20092 iHopkins et al.ll2010l ). Minor mergers are plausibly 
the driver of size growth, given that ma jor mergers are 
less e fficient in 'puffing up' galaxies fe.g.. lBezanson et"al"1 
2009), such that major merging would overproduce the 
numb er of massive galax ies in the present-day universe 
(e.g.. lMcLure et al~ll2012l ). Given the strong evidence for 
contin uous evolution of the number of ear ly- type galaxies 
(e.g., IBell et all 12003: iFaber et alJl200l and size evolu- 
tion, one may expect that their structural properties also 
evolve, especially if merging is invoked as an evolution- 
ary mechanism. To test this, H12 compared the axis- 
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ratio distribution of a large sample of early-type galaxies 
at z ~ 0.7 with that of the local population, but found 
no evidence for shape evolution for masses larger than 
3 x 10 10 Ai Q . However, there could be shape evolution 
at higher redshifts, as major merging occurred mor e fre- 
quently at early epochs (e.g., iRobaina et al.ll201Q[ ) and 
the early-type galaxy number was much lower at z > 1 
t han it is today- 
Ivan der Wei et all (|2011l hearafter vdWll) showed 
that a good portion of a small sample of z ~ 2 early- 
type galaxies appeare d flat in project i on, in dicative of 
a disk-like structure. iChevance et al.l (|2012f ) point out 
that existing, samples are too small to confirm or rule 
out evolution in the axis ratio distribution, but they at- 
test, based on the Sersic index distribution, that these 
galaxies have more disk-like structural pro perties than 
presen t-day early-type galaxies. Recently, iBruce et al.1 
(2012) used bulge-disk decompositions of massive early- 
type galaxies at z > 1 to show that man y of them host 
pronounced disks. IWhitakc r et al.l (|2012t ) show a hint of 
an evolving axis ratio distribution of early-type galaxies 
out to z ~ 1.5, but these authors did not explore this in 
detail. 

In this paper we explore the (projected) shapes for a 
large sample (^400 objects) of early-type galaxies, se- 
lected to have masses log(M/M & ) > 10.7 and low star 
formation rates (sSFR < l/3tii(z), see Section 2.5.), 
which were drawn from a wide, high-resolution, near- 
infrared (if -band) imaging mosaic from VLT/HAWK-I; 
we investigate whether early-type galaxies at z > 1 show 
evolution in structure compared to present-day counter- 
parts. Specifically, we will address the question whether 
z > 0.8 early-type galaxies are more or less disk-like than 
at the present epoch, and whether early-type galaxies at 
those epochs also become rounder with increasing mass, 
as seen today. 

The structure of this paper is as follows. In § 2 we de- 
scribe the data and select our sample of early-type galax- 
ies. In § 3 we analyze the projected axis-ratio distribution 
and its evolution since z ~ 2. In § 4 we summarize our 
conclusions. 

In this paper, we use AB magnitudes and adopt the 
cosmological parameters (fiM^A^)=(0-3,0.7,0.7). 

2. DATA 

The first step is to compile a catalog with photometric 
redshifts and stellar masses of galaxies in the extended 
Chandra Deep Field South (ECDFS) fr om MUSYC 
(Mul ti wavelength Survey by Yale- Chile, iTavlor et al.1 
120091 : iCardamone et all 120101 ). Then we use high- 
resolution VLT/HAWK-I K-band imaging available over 
essentially the full E-CDFS (ESO Program ID: 082.A- 
0890) to determine structural parameter (sizes, Sersic 
indices and projected axis ratios) for these galaxies. The 
high fidelity of these measurements, verified through the 
comparison with results from HST imaging, allows us to 
select early- type galaxies up to z — 1.8. 

2.1. Multi-Wavelength Data and SED Fitting 

MUSYC compiled observations in 32 bands, ranging 
from the UV to the ne a r-infr ared for the ECDFS, for 
which ICardamone et all (|2010l ) provide an optically se- 
lected catalog that we use here. We use the method 



and algorithms described bv lWuvts et al.l ()2011aD to in- 
fer photometric redshifts, stellar masses and rest-frame 
colors. Briefly, to estimate photometr ic redshifts (z p hoto) 
we use E Az Y (jBrammer et al.1 120081 ) , and to estimate 
stellar masses, star - formation rates, and rest-frame col- 
ors we use FAST (jKriek et al.l 120091 ). We only include 
objects with significant detection of J, H and if-band 
imaging, and reject stars by choosing only objects with 
J-K > 0.0 5 We adopt thelBruzual fc Chariot! ([2003f) 
model, and a Chabrier' (|2003f ) stellar Initial Mass Func- 
tion. A range of ages, star formation histories and extinc- 
tion parameters is explored. This parent catalog contains 
19642 objects. 

HI 2 independently determined the stellar masses of 
early-type galaxies in the ECDFS in the redshift range 

0. 6 < z < 0.8. Those mass estimates are designed to 
match the stellar mass estimates of present-day early- 
type galaxies. Since we aim to do the same we add 0.1 
dex to all our stellar mass estimates to correct for the 
median difference between the galaxies that are included 
both in our parent sample and the H12 sample. This cor- 
rection is likely incorrect for star-forming galaxies, but 
those are not considered in this work. A full investiga- 
tion of the absolute mass scale for z > 1 is beyond the 
scope of this paper. 

Star formation rates (SFRs) were derived following the 
procedures outlined in iWuvts et al.l (|2011bf ) . Briefly, the 
unobscured SFR traced by the UV was added to the 
dust- reemitted SFR inferre d from FIDEL 24 /im photom- 
etry (jMagnelli et al.1 120091 ) for 24 /mi-detected sources, 
and for sources without 24 /im detection a dust-corrected 
SFR was derived from stellar population modeling of the 
U-to-8 /im SED. 

2.2. High-Resolution, Near-Infrared VLT/HAWK-I 
imaging 

High- resolution if -band imaging from VLT/HAWK- 

1, is central to our study to provide the structural pa- 
rameters. We have obtained 1-hour exposures for each 
of 16 adjace nt tiles in a 30'x3 0' mosaic that covers the 
full ECDFS (jTavlor et al.ll2009l ). which is coincident with 
the HST/ ACS coverage from GEMS (Galaxy E volution 
from Morphologies and SEDs, iRix et al.ll200l . Obser- 
vation and reduction of the HAWK-I images have been 
performed by S. Z. using a customized pipeline based on 
the original version distributed by ESO. Most notably 
we implemented improved recipes for the construction of 
the master flat field and for the frame coaddition, which 
properly take into account object masks and variance 
maps. The new effective mask implementation in partic- 
ular eliminates the effects of background over-subtraction 
which is often seen in correspondence of (bright) sources. 
The 5-cr point source limit is if (AB)=24.3, and the Point 
Spread Function (PSF) has a FWHM of 0.5" or smaller 
across the field. At that resolution, this ground-based 
if -band imaging allows us to quantify the rest-frame op- 
tical structural properties of z > 1 galaxies with spa- 
tial resolution that differs by no more than a factor of 
~ 2 from that obtained with HST/WFC3 in the ii-band 
(see Figure [T]) . The image quality of our data is sub- 
stantially better than the ^0.6-0.7" seeing data used by 
IWhitaker etail (|20ll) . 

2.3. Galaxy Structural Parameters 
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Fig. 1. — Comparison betwee n VLT/HAWK-I imagin g; (top left) used for the analysis in this paper HST/WFC3 imaging (top right; 
Early Release Science data from Windhorst et al. 2011) of the same galaxies. The bottom panels show zoomed-in versions of the images 
of the HAWK-I and WFC3 images of 6 early-type galaxies in the redshift, stellar mass range of interest and circularized half-light radii of 
HAWK-I, as indicated. Each panel is about 2" on the side. These examples indicate that the resolution and depth of the HAWK-I imaging 
are sufficient to infer the projected axis ratio q, consistent with the values obtained from the WFC3 images. 
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TABLE 1 

Massive early- type galaxies: 0.6 < z < 1.8, 
loq(M/M Pl ) > 10.7 AND sSFR <1/Mh {z). ID IS 
THE SAME AsICARDAAIONE ~T AL.I l|20TOI) , REDSHIFT 
( z phot) AND STELLAR MASS (M) ARE FROM SED 

FITTING. K-BAND MAGNITUDE(K) , EFFECTIVE 
RADIUS (Re) AND PROJECTED AXIS RATIO (q P roj) 
ARE FROM THE GALFIT RESULTS. 



ID 


z phot 


M [logM Q ] 


K [mag] 


Re 


Qproj 


4032 


1.15 


10.72 


20.68 


0.25 


0.49 


4213 


1.11 


10.73 


20.25 


0.35 


0.83 


4619 


1.00 


10.87 


21.72 


0.09 


0.87 


4844 


1.59 


10.94 


20.94 


0.35 


0.56 


5375 


1.13 


11.33 


19.46 


0.44 


0.96 



We use GALAPAGOS (Galaxy Analysis over Large Ar- 
eas: Parameter Assessment by GA LFITting Objects 
from SExtractor, iBarden et al.ll2012h to separately pro- 
cess each of the 16 HAWK-I tiles. Here we briefly de- 
scribe the process as r elevant for the presen t study. For 
a full description, see IBarden et all (|2012f ). GALAPAGOS 
first constructs a c atalog with SExtractor (Software for 
source extraction, Bertin fc Arnoutsl 119961 ). We choose 
the SExtractor detection parameters such that the cat- 
alog is only complete down to K(AB)=23: as we will 
demonstrate below this is well beyond the limit down 
to which structural parameters can accurately be deter- 
mined. GALAPAGOS then creates image and noise cutouts 
for each object, including neighboring objects as neces- 
sary. The noise map is obtained from the variance map 
that was produced for each of the HAWK-I tiles. The 
background is estimated for each object by identifying a 
set of sky pixels that are not influenced by any of the 
objects in the catalog. 

For each of the 16 tiles a single star is taken as the 
PSF, chosen among the 5 brightest, isolated stars in each 
tile. This choice is made after subtracting a flux-scaled 
version of each of those 5 stars from ~25 stars in a tile 
and examining the residuals from this fit. The star that 
produces the cleanest residuals is selected as the PSF for 
that tile. 

Then GALFIT (v3.0.3, iPeng et alJl2lM 1201 Oaf) is called 
to perform the actual measurement of the structural pa- 
rameters. A single Sersic profile is fit to each target ob- 
ject. Neighboring objects are either masked or fit simul- 
taneously. The free parameters in the fit are position, 
magnitude {m), effective radius as measured along the 
major axis (R e ), Sersic index (n), axis ratio (q), and po- 
sition angle. The input values of these parameters are 
taken from the SExtractor catalog (with the exception 
of n, for which 2.5 is adopted). 

2.4. Precision and Accuracy of Axis-Ratio 
Measurements 

Our ability to determine axis ratios for distant galax- 
ies is the limiting factor in our study. To establish the 
precision and accuracy of our measurements we com- 
pare the axis ratios infe rred from the HAWK -I imaging 
with those from GEMS (|Haussler et al.|[2007t ) for galax- 
ies in the redshift range 0.6 < z < 0.8 (the sample 
from H12). The wavelength difference between GEMS 
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Fig. 2. — Comparison of axis-ratio measurements from HAWK-I 
/f-band and from HST imaging. The crosses represent the compar- 
ison with measurements from WFC3 F160W imaging for a HAWK- 
I K-b&nd selected sample (regardless of galaxy type and redshift). 
The plus signs represent the comparison with measurements from 
HST/ACS F850LP of the H12 sample, which includes early-type 
galaxies in the redshift range 0.6 < z < 0.8. The top panel di- 
rectly compares the VLT and HST axis-ratio measurements; the 
bottom panel shows the difference between the VLT and HST mea- 
surements as a function of HAWK-I _ft"-band magnitude. The red 
and orange bars represent the median and standard deviation for 
a series of magnitude bins. The standard deviation represents the 
measurement uncertainty in the VLT-inferred axis ratio, assuming 
the HST-based value as 'truth'. For galaxies, regardless of type, 
brighter than K = 22 the uncertainty is smaller than or equal to 
0.1, and we adopt this as the magnitude limit for our study. Inter- 
estingly, the accuracy of the axis-ratio measurement is good down 
to at least K ~ 24: systematic effects in the axis-ratio measure- 
ments are small compared to the random uncertainty. 



(z-band) and HAWK-I (.fT-band) could cause intrinsic 
differences between axis-ratio measurements that cannot 
be attributed to measurement errors. For this reason, we 
also compare with the axis ratio estimates from objects in 
our sample that are contained with HST/W FC3 F160W 
imaging from ERS (jWindhorst et all 120111) . GALAPAGOS 
is deployed in a similar fashion as described above - 
full details will be provided by van der Wei (2012, in 
prep.). Through the comparison with HST/ACS and 
HST/WFC3 (see Figured]) we find that the precision of 
our HAWK-I axis ratio estimates is better than 10% for 
galaxies brighter than K(AB)=22. In addition, we see 
no systematic difference between HST and VLT mea- 
surements shown in Figure [2] The differences in the me- 
dian are -0.019 and -0.025, respectively. Our axis-ratio 
measurements remain accurate over the entire magnitude 
range of our sample. In Figure [3] we show that, in addi- 
tion, the accuracy of our axis ratio measurements does 
not depend on redshift. 
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Fig. 3. — Comparison of axis-ratio measurements from 
VLT/HAWK-I and HST/WFC3 imaging as a function of redshift. 
The difference in axis ratio measurements, qu AWK-I —1WFC3 has 
been multiplied by the total K band flux in units of the K band 
flux for a K = 21 object, such that variations in magnitude and 
their effect on the measurement uncertainty are accounted for. For 
galaxies in a narrow magnitude range (here, 20 < K < 22), there 
is no systematic redshift dependence: the uncertainty for galaxies 
at 0.5 < z < 1.0 is comparable to the uncertainty for galaxies at 
z > 1. The trends shown in Figure [2] and this figure imply that 
the accuracy of our axis ratio measurements is 10% of better for 
all galaxies in our sample. 

2.5. Sample Selection 

We match the MUSYC-based catalog (containing pho- 
tometric redshifts and stellar population properties) with 
the HAWK-I based catalog (containing iv-band mag- 
nitudes and structural parameters) by searching within 
apertures with radius 1" . In Figure S] we show the dis- 
tribution of if-band magnitude as a function of stellar 
mass in three redshift bins. Given the magnitude limit 
of K = 22 that we adopted to ensure precise axis-ratio 
measurements (see above), we find that our catalog is 
complete down to log(M/M®) — 10.7 for all redshifts 
z<1.8. 

We select galaxies with stellar masses log(M/M@) > 
10.7 and redshifts 0.6 < z < 1.8, using spectroscopic red- 
shifts when available, and, otherwise, photometric red- 
shifts. We select as early-type galaxies those with spe- 
cific star- formation rates (sSFR — SFR/Mq) smaller 
than l/3tfr(z), where tji(z) is the Hubble time at the 
photometric redshift of the galaxy. This is similar to the 
strategy adopted by HI 2, who also selected their samples 
by star-formation activity, rather than morphological ap- 
pearance. 

Visual morphological classifications are not possible 
for the z > 1 galaxies in our sample as they are too 
faint. Automated classifiers based on, for example, con- 
centration or Sersic index have the problem that, by def- 
inition, they will select against disk-like objects. This 
motivates our choice to select early-type galaxies by 
their star-formation activity. This is justified further 
by the well-established agreement between morphologi- 
cal appeara nce and star form ation activity at all redshifts 
z < 1 (e.g., iBell et a l. 2004), and the correlation between 
structure and star formation activi ty that is observed to 
persist out to at least z ~ 2 (e.g., iSzomoru et alJl2011t 
IPatel et al]l20H . 

In Figure [S] we show that our selection technique 
is compatible with the rest-frame UVJ selection tech- 
nique that is o ften adopted to identify passive, early- typ e 
galaxies (e.g, W uvts et al.ll2007t iWilliams et all 120090 . 
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Fig. 4. — HAWK-I X-band magnitude versus stellar mass in three 
redshift bins. Early- type galaxies, selected by their low sSFR (see 
Section 2.5.), are indicated by red symbols. The limiting factor 
in selecting our sample is set by the axis ratio estimate precision 
(see Figure 0, which sets our magnitude limit to K = 22 (the 
horizontal lines). The implication is that in our highest redshift 
bin we are complete down to a mass limit of log(M / 'Mq) ~ 10.7 
(the vertical lines). 



Almost all of our low sSFR galaxies would also be identi- 
fied as early-type galaxies by selecting them in this UVJ 
diagram. There is a substantial number of galaxies in the 
'passive' UVJ color-color box that are star forming ac- 
cording to our direct star formation estimates fits. This 
number is far larger than the number of galaxies with low 
star formation rates that are located outside the 'passive' 
UVJ color box. While we include the latter in our sub- 
sequent analysis, we note that this does not affect our 
results. Samples selected by UVJ, star-formation rate 
or a combination all have have the same median axis 
ratio within 0.05. 

In Figure[S]we show the correlation between rest-frame 
U — V color and sSFR. Even in our highest redshift bin 
the populations of star-forming and passive galaxies sep- 
arate cleanly, which implies little cross-contamination be- 
tween the two types of galaxies. 

Finally, we visually inspected all GALFIT fitting results, 
rejecting 10 objects with corrupted fits. The final sam- 
ple consists of 394 galaxies (see Table ^ in the red- 
shift range 0.6 < z < 1.8, log(M/M Q ) >= 10.7, and 
sSFR < l/3t H (z). There are 134, 163 and 97 galaxies 
in the redshift bins, 0.6 < z < 0.8, 0.8 < z < 1.3 and 
1.3 < z < 1.8, respectively. 
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TABLE 2 
Statistical properties of q pr0 j 



log(M/M e ) > 10.7 


number 


median 


mean 


stddev 


median of z 


1.53 < z < 2.31 (vdWll) 


14 


0.665 


0.681 


0.186 


1.69 


1.3 < z < 1.8 


97 


0.634 


0.612 


0.200 


1.38 


0.8 < z < 1.3 


163 


0.646 


0.619 


0.213 


1.09 


0.6 < z < 0.8 


134 


0.734 


0.662 


0.222 


0.69 


0.6 < z < 0.8 (H12) 


533 


0.710 


0.692 


0.179 


0.69 


0.04 < z < 0.08 (vdW09) 
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0.181 


0.06 
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Fig. 5. — Rest-frame U-V color versus rest-frame V-J color, a di- 
agnostic diagram to distinguish passive and star-forming galaxies. 
The early-type galaxies in our sample, selected by having SED- 
based sSFR(=SFR/M)< l/3t£r(z), are indicated by red symbols. 
Th ese are essentia lly always located in the 'passive box' as defined 
by Williams et al. (2009)as indicated by the polygons. That our 
selection criteria are conservative is indicated by the presence of a 
substantial number of objects with U-V and V-J colors consistent 
with those of passive galaxies but with sSFR larger than our limit. 

3. STRUCTURAL EVOLUTION OF MASSIVE EARLY-TYPE 
GALAXIES 

We now examine two aspects of the structural evolu- 
tion of early-type galaxies. First, we address the question 
whether the most massive galaxies are intrinsically round 
at redshifts 1 < z < 2, as is observed at z < 1 by vdW09 
and H12. Second, we address the question whether, in 
general, early- type galaxies are flatter or rounder (e.g., 
disk- like or bulge-like) at z > 1 than at z < 1; H12 
showed that there is little evolution up to z < 1, but 
vdWll presented tentative evidence for a higher inci- 
dence of disk- like early- type galaxies at z ~ 2. 



Fig. 6. — The correlation between rest-frame U — V color and 
sSFR for massive galaxies in the three redshift bins. The red sym- 
bols represent galaxies in the 'passive' UVJ color-color box (see 
Figure 0; the black symbols represent star- forming galaxies. Our 
selection criterion for early-type galaxies is based on a sSFR cut 
(sSFR < l/3f, j/(z)), here illustrated by the horizontal, dotted lines 
which correspond to the sSFR cuts in the center of the redshift bins. 
For all redshift bins the two types of galaxies are cleanly separated; 
moving the sSFR cut by a modest amount does not strongly affect 
our selection of early-type galaxies. Note that aliasing occurs for 
low-sSFR objects as a result of template choices in the SED fitting. 

3.1. The Mass-Dependence of Early- Type Galaxy 
Shapes up to z ~ 2 

In Figure [7] we show the projected axis ratios of our 
sample of early-type galaxies as a function of their stel- 
lar mass, split into three redshift bins, each with ~ 100 
galaxies. 

Qualitatively speaking, Figure [7] shows the same trend 
in all redshift bins: there appears to be a mass depen- 
dence for the projected axis-ratio, such as was seen in 
vdW09. As previously demonstrated by H12, at z ~ 0.7 
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Fig. 7. — The projected axis ratio versus stellar mass for our 
early- type galaxy sample, split into three redshift bins. The sym- 
bols arc the HAWK-I ii"-band images of the galaxies. The lines rep- 
resent running median values of the axis ratio. Up to our highest 
redshift bin, the most massive galaxies appear to be the roundest 
(see also Figure[8]and Sec 3.1). Overall, the galaxies in the highest 
redshift bins are also flatter than their lower-redshift counterparts 
(see also Figure|9]and Sec 3.2). 



we find that the most massive (log(M/M®) > 11.3) 
early-type galaxies are round in projection and we find 
a lack of objects that are flat in projection among the 
most massive early- type galaxies. In Figure [7] a simi- 
lar trend is seen for early- type galaxies at z > 1: at 
and below masses of log{M /Mq) ~ 11 the galaxies in 
our sample show a broad range in axis ratios, whereas 
more massive galaxies are predominantly round. This is 
consis tent with the conclusion reached bv lTargett et al.l 
(|2011D . who found that very luminous radio galaxies at 
z ~ 2 have morphological properties similar to today's 
most massive elliptical galaxies in clusters. 

In Figure |5] we compare the axis-ratio distributions of 
galaxies with masses below and above log(M /Mq) ~ 
11.3. The Mann- Whitney U test shows that the trend 
seen in Figure [7] is marginally significant in each of the 
redshift bins. For the combined sample of galaxies at 
rcdshifts 0.8 < z < 1.8 we find that the most massive 
galaxies are rounder than the less massive galaxies with 
high confidence (P = 2.51 x 10~ 4 , i.e., a 3.66cr result). 
For comparison, the equivalent value from the classifical 
Kolmogorov-Smirnovtest is also low (Pks — 1-81 X 10~ 3 , 
i.e., a 3.12a result). 

Combining the results from vdW09, H12 and from 
this paper, we conclude that at all redshifts < z < 2 
the most massive early-type galaxies are predominantly 
round in projection. This implies that at all these cos- 
mological epochs the formation mechanism for such mas- 
sive galaxies precludes the formation of or requires the 
destruction of pre-existing disks. Our interpretation is 
that merging, accompanied with little dissipation and 
star formation, has been the dominant formation chan- 
nel for galaxies more massive than log(M /Mq) = 11.3 
since z ~ 2. This process is rep roduced in early- 
type galaxy for mation models (e.g., iNaab et al.l 120091 : 
IQser et al.ll20Tl . and the observed major merger rates 
up to z ~ 2 are consis tent with theoretical expectations 
(e.g- IMan et~aTl[2012h . 

3.2. Shape Evolution at 1 < z < 2 

Now we turn to the question whether the structural 
population properties evolve with redshift. vdWll found 
an indication of a high incidence of disk-like early-type 
galaxies at z ~ 2, but their sample is too small to confirm 
or rule out evolution with respect to the structure of 
early-type galaxies at the present day. In Figure [9] we 
compare the shape distributions of early-type galaxies 
with masses log(M /Mq) > 10.7 across a broad range in 
redshift, < z < 1.8. We reproduce the lack of strong 
e volution in the range < z < 0.8, reported by H12. 

iVulcani et al.l ([201 ll ) showed that cluster early-type 
galaxies are rounder at z=0.5-l than in the local uni- 
verse. This result is not necessarily at odds with our mea- 
surements: their signal is mostly driven by galaxies below 
our mass limit and, moreover, struct ural differences be- 
tween cluster and field galaxies (e.g., Ivan der Wei et all 
12011) can be explained by environmental processes such 
gas stripping that can produce flat early-type galaxies 
in clusters. An indication that the chosen mass range is 
relevant is that low-mass galaxies show stronger environ- 
mental dependenc ies in their proper ties than high-mass 
galaxies (see, e.g. JPeng et al.ll2010b[ ). Observations over 
a broader range of environments and down to lower stel- 
lar masses than what is possible with our data set are 
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Fig. 8. — Axis ratio histograms of early-type galaxies with masses 
10.7 < log(M I Mq) < 11.3 (dashed lines) and early-type galaxies 
with masses log(M / Mq) > 11.3 (solid lines). The bottom panels 
show the distributions (df /dq pro j) of different redshift bins, as in- 
dicated; the top panel combine all redshift bins and show the cumu- 
lative histograms as a function of projected axis ratio (/(< q pr oj), 
where the color coding corresponds to redshift, following the color 
coding in the bottom panels. The top panel also shows in black 
the SDSS-based sample of local early-type galaxies from H12 in 
the same mass range. At all redshifts the most massive galaxies 
galaxies are the roundest; the significance of this observation is 
confirmed by the M-W statistical test (as indicated by the listed 
probabilities that the samples are statistically the same). 

required to address these issues. 

However, in our higher redshift bins (z > 0.8) we see 
an excess of fiat galaxies and a detectable evolution com- 
pared to the z < 0.8 samples. A Kolmogorov-Smirnov 
(K-S) and Mann- Whitney U (M-W) tests shows that this 
trend is highly significant (see Tables [3] and U]). Note that 
the axis ratio distribution from vdWll, due to their small 
sample size, is consistent with both the low- and high-z 
axis ratio distributions we analyze here. 

Figure [TU] summarizes the observed evolution in the 
mean projected shape of early type galaxies. Beyond 
z~lwe see a gradual decrease in the median axis ratio 
of early-type galaxies, which is the result of an increased 
fraction of galaxies with axis ratios smaller than 0.6. The 
observed axis-ratio distributions imply that the intrinsic 



Fig. 9. — Axis ratio histograms (df /dq pro j) of early-type galaxies 
with all masses log(M / Mq) > 10.7 (our full sample) as a func- 
tion of redshift in bottom three panels, always compared with the 
present-day sample from H12. The top panel shows cumulative 
histograms for all redshift bins as a function of projected axis ra- 
tio /(< q pro j). Up to z=0.8, as confirmed by the M-W statistical 
test (of which the probability is given that the samples are statis- 
tically the same) there is no significant redshift evolution in the 
projected axis-ratio distribution, consistent with the results from 
H12. At z > 1, we find that early-type galaxies are flatter than 
their present-day counterparts. 

thickness of the typical early-type galaxy is no more than 
about 0.4. In a forthcoming paper we will quantify this 
intrinsic shape through detailed modeling. 

The observations presented here can now put on a 
firmer footing the claim by vdWll that many high- 
redshift early-type galaxies display a disk-like struc- 
ture and that, plausibly, these galaxies will grow into 
larger, intrinsically round galaxies through merging. The 
vdWll sample had been too small to confirm or rule 
evolution at a fixed stellar mass, which is important in 
order to differentiate evolutionary paths that galaxies of 
a given mass take at different cosmic times. Even at the 
present-day, the typical early-type galaxy with a mass be- 
low log(M/M Q ) = 11 is rather flattened (vdW09), which 
also corresponds to a disk-like ki nematic structure (rota - 
tion) as shown by, for example, Emscllc m et al.l (|2011l) . 
Here we show that such galaxies had yet thinner intrinsic 
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Fig. 10. — Top panel: median axis ratio of early-type galaxies more massive than log(M / Mq) > 10.7 as a function of redshift. The 
sample presented in this paper is represented by the colored data points; the black data points are taken from the literature as indicated. 
The error bars are estimated by bootstrapping. Bottom panel: fraction of galaxies with axis ratios smaller than indicated, again as a 
function of redshift. Beyond z ~ 1, early- type galaxies typically have smaller axis ratios than present-day counterparts in the same mass 
range, indicative of a more disk-like structure. All error bars are inferred by 95% confidence level compare to the same fraction. 



TABLE 3 

KOLMOGOROV-SmIRNOV TEST: SIGNIFICANT LEVEL (log(M /Mq) > 10.7) 



log(M/M e ) > 10.7 


1.53 < z < 2.31 


1.3 < z < 1.8 


0.8 < 2 < 1.3 


0.6 < z < 0.8 


0.6 < z < 0.8 


0.04 < z < 0.08 


1.53 < z < 2.31 (vdWll) 




0.463 


0.549 


0.671 


0.668 


0.688 


1.3 <z < 1.8 


0.463 




0.875 


*1.02E-2*ta] 


*2.58E-3* 


*1.54E-4* 


0.8 < z < 1.3 


0.549 


0.875 




♦1.26E-2* 


*8.93E-4* 


*1.21E-4* 


0.6 < z < 0.8 


0.671 


*1.02E-2* 


*1.26E-2* 




0.060 


0.035 


0.6 < z < 0.8 (H12) 


0.668 


*2.58E-3* 


*8.93E-4* 


0.060 




0.584 


0.04 < z < 0.08 (vdW09) 


0.688 


*1.54E-3* 


*1.21E-4* 


0.035 


0.584 





a Star symbol (*) represent the significant probability is smaller than 5%. It implies the distributions arc distinguishable. 



TABLE 4 

Mann- Whitney test: significant level (log(M/MQ) > 10.7) 



log{M/M & ) > 10.7 


1.53 < 2 < 2.31 


1.3 < 2 < 1.8 


0.8 < 2 < 1.3 


0.6 < 2 < 0.8 


0.6 < 2 < 0.8 


0.04 < 2 < 0.08 


1.53 < z < 2.31 (vdWll) 




0.139 


0.167 


0.456 


0.372 


0.373 


1.3 < z < 1.8 


0.139 




0.315 


*1.12E-2* la] 


*1.20E-4* 


*3.00E-5* 


0.8 < z < 1.3 


0.167 


0.315 




*2.06E-2* 


*8.71E-5* 


*9.18E-6* 


0.6 < 2 < 0.8 


0.456 


*1.12E-2* 


*2.06E-2* 




0.249 


0.234 


0.6 < 2 < 0.8 (H12) 


0.372 


H.20E-4* 


*8.71E-5* 


0.249 




0.445 


0.04 < z < 0.08 (vdW09) 


0.381 


*3.00E-5* 


*9.18E-6* 


0.234 


0.445 





a Star symbol (*) represent the significant probability is smaller than 5%. It implies the distributions arc distinguishable. 
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shapes at z > 1. Good correspondence between rotation 
of the stellar body and flattening has been confirm ed up 
to z ~ 1 by Ivan der Wei k van der Marell (|2008D . Al- 
though such confirmation is, so far, lacking at higher 
redshifts it is reasonable to assume that flattening im- 
plies rotation at all redshifts. 

We note that beyond the simple observation that early- 
type galaxies, here defined as galaxies with little star- 
formation activity, are on average flatter at z > 1 than at 
lower redshifts, we cannot distinguish between the differ- 
ent varieties of such galaxies - such as quiescent spirals, 
barred/ringed SOs, etc. - and their evolution. 

Our results do not argue against a morphological tran- 
sition that occurs along with the t runcation of star for- 
mation, as argued by, for example, iBell et all (|2012f) on 
the basis of the high Sersic indices of quiescent galaxies: 
high-rcdshift early-type galaxies such as those studied 
here, despite the generally disk-like character inferred 
from their axis ratio distribution, are not pure disks 
and have higher Sersic indice s than star forming galax- 
ies (e.g., IWuyts et al.l l2011af). B ased on bulge-disk de- 
compositions. lBruce et al.l(|2012f) arrive at the conclusion 
that the early-type galaxy population at z > 1 is mix of 
bulge- and disk-dominated galaxies, indicating that the 
transition from actively star-forming to quiescent need 
not always coincide with the formation of a dominant 
bulge. Taken together, these observations are consistent 
with a picture in which gas had time to settle in a disk 
before star formation was truncated, and that this disk 
wholly or partially survived the process that truncated 
star formation. 



4. CONCLUSIONS 

We measured the projected axis ratios from 
VLT/HAWK-I if-band imaging of a sample of early- 
type galaxies in the redshift range 0.6 < z < 1.8 selected 
by their low specific star formation rates. We find 
that at all redshifts z < 2 the most massive galaxies 
log(M /Mq) > 11.3 are predominantly round. The 
lack of very massive, highly flattened galaxies suggests 
a universal ceiling mass for the formation of disks, 
independent of cosmic epoch. In order for galaxies to 
grow beyond this ceiling mass, separate evolutionary 
channel, presumably merging, has to be invoked. 

In the full sample (log{M/M Q ) > 10.7), we find, at all 
redshifts, a large range in projected axis ratios, reflecting 
a more disk-like structure. In addition, we find quantita- 
tive evidence that early-type galaxies at z > 1 are more 
disk-like than their equally massive, present-day coun- 
terparts. Therefore, for most early-type galaxies gas had 
time to settle into a disk before star formation ceased to 
produce a more passive galaxy. Plausibly, such galaxies 
grow in mass over time, through mostly dissipationless 
merging and accretion of satellites, losing some of the 
angular momentum, and growing in size. 



Y.-Y. Chang was funded by the IMPRS for Astronomy 
& Cosmic Physics at the University of Heidelberg and 
the Marie Curie Initial Training Network ELIXIR of the 
European Commission under contract PITN-GA-2008- 
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